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ABSTRACT 

We show that the removal of angular momentum is possible in the presence of large scale magnetic 
stresses in geometrically thick, advective, sub-Keplerian accretion flows around black holes in steady- 
state, in the complete absence of a-viscosity. The efficiency of such an angular momentum transfer 
could be equivalent to that of a-viscosity with a = 0.01 — 0.08. Nevertheless, required field is well 
below its equipartition value, leading to a magnetically stable disk flow. This is essentially important 
in order to describe the hard spectral state of the sources, when the flow is non/sub-Keplerian. We 
show in our simpler 1.5-dimensional, vertically averaged disk model that larger the vertical-gradient of 
azimuthal component of magnetic field, stronger the rate of angular momentum transfer is, which in 
turn may lead to a faster rate of outflowing matter. Finding efficient angular momentum transfer, in 
black hole disks, via magnetic stresses alone is very interesting, when the generic origin of a-viscosity 
is still being explored. 

Subject headings: accretion, accretion disks — MHD — jets and outflows — X-rays: binaries — 
galaxies: active 


1. INTRODUCTION 

iBlandford fc Pavne I (|1982ll showed the possibility of energy and angular momentum removal from a Keplerian 
accretion disk by the magnetic field lines that extend from the disk surface to large distance. In the framework of infinite 
conductivity and self-similar model, they showed that the magnetic stresses can extract the angular momentum from a 
geometrically thin accretion disk, which helps matter to accrete, independent of the presence of viscosity. Furthermore, 
they argued that such a mechanism is responsible for the observed jets/outflows from accreting sources, when magnetic 
stresses convert a centrifugal outflow into a collimated jet. The disk matter has been argued to be outflowing through 
the outgoing field lines. The time evolution of axisymmetric, weak magnetic fields threading geometrically thin, 
Keplerian acc r etion disks with finite conductivity in a specific model framework was furthermore investigated by 
iLubow et al. I (Il994l l. however without considering possible angular momentum transfer by the magnetic field. On 
the other hand, in the presence of infinite conductivity, the magnetic field, in the same model framework which does 
not consider the contributio ns from the m agnetic st resses, would be kept on amplifying b y the accretion of gas , till it 
stops the accretion (also see lSDruit1l20I3ll . However. iBisnovatvi-Kogan fc Lovelace I (l2007ll and lLovelace et akl (l2009f) 
showed that the radially inward flow is possible for plasma- /3 > 1 and Prandtl number > I, in the stationary channel- 
type flows having small optical depth in the absence of turbulent viscosity, which also could exhibit electromagnetic 
outflows for smaller Prandtl number. They showed that the large-scale field keeps drifting inward until a stationary 
state arises, when the magnetic, centrifugal, and gravitational forces become comparable. This furthermore reveals 
the flow velocity profile differing significantly from the Keplerian profile. 

The idea of exploring magnetic stress in order to explain astrophysical systems was, in fact, implemented much 
earlier. For example, the solar wind was understood to have decreased Sun’s angular momentum through the effect 
of magnetic stres ses (see, e,g.. IWeber fc Day s I Il967ll . the proto-stellar gas clouds m ight have been co n tracte d by 
magnetic effects (iMouschovias fc Paleologou 11198011 . In the context of accretion disk, lOzernov fc Usov I (Il973ll and 
IBlandford (|I976ll showed that t he energy is possible to be extracted continuously by electromagnetic torques and 
twisted field lines. Furthermore, ICao fc Spruit I (|2002t) showed, by a linear stability analysis of the accretion disks, 
that angular momentum is possible to be removed by the magnetic torque exerted by a centrifugally driven wind. The 
same authors (|Cao fc Spruit ll2M^ also discussed that moderately weak fields can cause sufficient angular momentum 
loss, via a magnetic wind to balance outward diffusion in geometrically thin accretion disks. However, plasma-/? has 
to be much smaller than unity to explain the t endency of strong flux bundles at the centers of disk to stay confined, 
as seen in numerical simulations. Nevertheless, iQgilvie fc Livio I (jI998ll showed, by solving the local vertical structure 
of a thin accretion disk threaded by a poloidal magnetic field, the shortcoming of launching an outflow and suggested 
for an existence of additional source of energy for successful launching of the outflow. 

However, observat ionally outflows/jets are mostly found to be emanating from the disk when it is in a hard state (e.g. 
iBelloni et al. ll^OOt) . which is non/sub-Keplerian. Note that jets appear to be highly heterogeneous with yelocities 
ranging from a few tens of million cm/sec to the escape yelocit y from the disk. Superluminal sources, however, appear 
to exhibit jet velocity around the speed of light (|Milev 11198011 . The jets are found in disks around stellar mass black 
hole sources (e.g. GRS I9I5-I-I05) as well as supermassive black hole sources (e.g. M87). 

Therefore, most of the modern models describing outflows/jets from the accretion disks are based on sub-Keplerian, 
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advective model when the flow has a significant radial velocit y, unlike the K e pleria n disk. For example, a class 


'.e pleria n i 

] IT^ . i 


in order to describe bipolar 


Jf:- 

of self-similar, advection dominated solutions was proposed by iNaravan fc Yi I 
outflows from black hole sources, e.g. Sgr A*. Later on, such a class of advection dominated s olution exhibiting 
outflo ws/jets was applied in many other contexts, e.g. cor e-collapsing disks and g amma-ray bursts (|Di Matteo et al. 1 
1200211 . low-radiative-efSciencv nucl ei of elliptical galaxies (|Di Matteo et al. 1 120021). 

In a different model framework, IChakrabarti fc TitarchulTI ( 19951 1. IChakrabarti I (|1999[ 1. IChakrabarti fc Manickarnl 
(|2000ll described advective, sub-Keplerian accretion flows in or der to explain outflows , quasi-periodic oscillations 
QPO s) and spectral states in black hole sources. Furthermore, iMukhopadhvav I (|2003ll . iMukhopadhvav fc Ghosh I 
200311 . iRaiesh &: Mukhopadhvav I (|2010all described general advective accretion flows (GAAFs) around black holes 
and neutron stars and showed the effects of rotation of the black hole on to the solutions. The last authors also 
included various cooling effects explicitly and showed how the solutions get affected by the cooling properties. 

However, all the abo ve models were formulated in the framework of Shakura-Sunyaev a-viscosity 
(jShakura fc Sunvaev I Il97,3l l. when the flows are assumed to be embedded with the plasma-/? >> 1. Hence, the 
matter tr ansport is assumed to be su pported via tur bulent viscosity, not by lar ge scale magnetic field, unlike that 
chosen bv iBlandford fc Pavne I (|1982fl . Nevertheless, iBhattacharva et al. I (l2010ll showed that the transport is also 
possible in the presence of outflow in a 2.5-dimensional accreting system; it does not matter whether the outflow is 
magnetic or hydrodynamic. Note that outflows and even jets can also be formed in the absence of magnetic field. 
This is likely to occur when the flow is radia t ion trapped and the accretion rate is super-Ed dington or super-critical 
(jLovelace et al. I[i994 iBegelman et al. II 2 OO 6 I: iFebrika Il2004 iGh osh fc Mukhopadhvav 11200^ . 


The 2.5-dimensional accretion model, proposed by 


Bhattacharva et al. I ( 201(1( 1. will be complete if the effects of (large 


scale) magnetic field is included therein. In that case, on e presumably can exp l ain ou tflow of matter plunging through 
the outgoing magnetic field lines more spontaneously, as IBlandford fc Pavne I (|1982ll did in the Keplerian framework. 
To the best of our knowledge, so far there is no attempt to obtain a self-consistent set of advective disk-outflow coupled 
accretion solutions in the presence of large scale magnetic field, which has lots of implication to explain low/hard state 
of sources. This, howeve r, has been discussed in some extent for circumstellar disks around young stars (see, e.g., 
iKbnigl fc Salmeron ll20Hll . without discussing the detailed solutions of all the dynamical variables. The present work 
steps forward in order to obtain such a set of solutions for black holes. 

In various numerical set ups, magnetohydrodynamic (MHD) simulations of accretion on to magnetized compact 
objects have alrea dy been explored. As ex amples, some of them considered axisymmetric systems in the presence 
of magnetosphere ((Romanova et aL~l 20111). some othe r aimed at investigating advection of matter and magnetic 
held in the turbulent/diffusive disks (|Dvda et al. II2013I) . when the held strength decreases due to reconnection and 
annihilation at a later time. Other groups, explored general relativistic magnetohydrodynamic (GRMHD) simulations 
of magnetically arrested accretion hows an d outhows around black holes, for toroidally and poloidally dominated 
magnetic helds (jTchekhovskov et al. l[2?)Tll ). They furthermore demonstr ated the possible extrac tion of net energy 
from a spinning black hole via the Penrose-Blandford-Znajek mechanism (iMcKinnev et al.~ll2012l) . Moreover, there 
were radiation MHD/GRMHD simulations of accretion and outhows around black holes, exploring three distinct 
how phases including the radiatively inefficient pha se which is similar to the hows considered in the present work 
(jOhsuga fc Mineshige II2011I: IMcKinnev et al. Il20l'4l . Some of the MHD simulations investigated the reasons behind 
the variability in low angul ar momentum, underluminous accretion hows in the vicinity of a supermassive black hole 
(jMoscibrodzka et al. 11200^ . However, all these works, to the best of our knowledge, considered the cases when any 
viscosity to be arisen from magnetorotational instability (MRI) leading effectively to the Shakura-Sunyaev a viscosity 
((Balbus fc Hawl^lI99Ill . 

Here we plan to investigate, semi-analytically, the effects of large scale magnetic held, with plasma-/? > 1 yet, on to 
the advective accretion hows in order to transport matter, however restricted to the simpler 1.5-dimensions. Therefore, 
we consider the how variables, averaged over the vertical coordinate, to depend on the radial coordinate only. While 
the vertical equilibrium assumption corresponds to no vertical component of velocity, we choose the vertical component 
of magnetic held to be non-zero. Although, in reality, a non-zero vertical magnetic held induces a vertical motion, 
in the platform of the present assumption, any vertical motion will be featured as an outward motion. Nevertheless, 
whether it is a vertical or outward transport, our aim is to furnish removal of angular momentum from the how via 
magnetic stresses, leading to the infall of matter towards black holes. 

The plan of the paper is the following. In the next section, we describe the set of magnetohydrody- 
namic/hydromagnetic equations, at the limit of very large Reynolds number, as is the case in accretion disks, de¬ 
scribing how model. This is basically the set of Navier-Stokes equation, but in the presence of magnetic shearing 
stresses (and Lorentz force), magnetic induction equation, the condition for the absence of magnetic monopole, and 
hnally the conservation of mass. Subsequently, we discuss the numerical solutions of the set of equations in §3 and its 
implications. Finally, we summarize the results along with a discussion in §4. 


2. MODEL HYDROMAGNETIC EQUATIONS 

We describe optically thin, magnetized, viscous, axi symmetric, advec t ive, v ertically averaged, steady-state accretion 
how, in the pseudo-Newtonian framework with the iMukhonadhyav I ()2002l) potential. The choice of the pseudo- 
Newtonian framework, for the present purpose, does not hinder any physics, compared to that would appear in the 
full general relativistic framework. Hence, the equation of continuity, vertically averaged hydromagnetic equations for 
energy-momentum balance in different directions are given by 

M = 4TTxphi9, (I) 
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( 5 ) 


when we assume that the variables do not vary significantly in the vertical direction such that d/dz —>■ Sijz ^ Si/h, 
when i = 1, 2,3, which is indeed true in the disk flows. Note that si, S 2 and S 3 are the degrees of scaling for the radial, 
azimuthal and vertical components of the magnetic field respectively. As a consequence, the vertical component of 
velocity is zero. Here M is the conserved mass accretion rate and the corresponding equation is the integrated version 
of the continuity equation, p is the mass density of the flow, d the radial velocity, P the total pressure including 
the magnetic contribution, F the force corresponding to the pseudo-Newtonian potential for rotating black holes, 
A the angular momentum per unit mass, Wxp the viscous shearing st ress written following the Shak ura-Sunyaev 
(|Shakura fc Sunvaev 111971^ prescription with appropriate modification (jMukhopadhvav fc Ghosh Il2003ll . h ^ z, the 
half-thickness of the disk, s the entropy per unit volume, T the (ion) temperature of the flow, and Q~ are the net 
rates of energy released and radiated out per unit volume in/from the flow respectively (when Qmagj Qvis^ Qmag 
are the respective contributions from viscous and magnetic parts). We furthermore assume, for the present purpose, 
the heat radiated out proportional to the released rate with the proportionality constants (1 — /^is) and (f — fm), 
respectively, for viscous and magnetic parts of the radiations. ri,r 3 indicate the p olytropic indices depending on 
the gas and radiation content in the flow (see, e.g., iRaiesh fc Mukhopadhvav Il2010al for exact expressions) and B^, 
B^ and Bz are the components of magnetic field. Note that, the independent variables x and z are the radial and 
vertical coordinates, respectively, of the flow, expressed in the units of GM/c^, where G is the Newton’s gravitation 
constant, M the mass of the black hole and c the speed of light. Accordingly, all the above variables are made 
dimensionl ess, e.g. d is expressed in the u nits of c. For any other details, e.g. model for see the existing 

literature (IRaiesh fc Mukhopadhvav Il 2010 allbl l. when 


Qvis Qvis 
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( 6 ) 


We furthermore do not consid er _ ^e heat generated and absorbed due to the nuclear reactions 

(jMukhopadhvav fc Chakrabarti 1 120001120011) . This is to emphasize that all the variables appearing in the equations 
are assumed to be their respective vertically averaged quantities. 

Hydromagnetic flow equations must be supplemented by (for the present purpose, steady-state) equations of induc¬ 
tion and no magnetic monopole, given by 


X/ X V X B + Vrn^^B = 0, 


( 7 ) 


d , . Bz , , 

— (xBx) + S 3 —= 0, (8) 

dx n 

when V and B are respectively the velocity and magnetic field vectors and is the magnetic diffusivity. On taking the 
ratio of the orders of the first to the second (diffusive) terms in equation ([7|), we obtain L\v\/vyn = Rm, when L being 
the order of the length scale of the system. Hence, when the Reynolds number, Rm, is very large, which is the case for 
accretion disks, the second term (which is associated with the magnetic diffusivity) in equation (jT]) can be neglected. 
However, this term can be rather important inside some localized regions in certain astrophysical systems due to subtle 
reasons. Nevertheless, for the present purpose, for simplicity we will neglect this term throughout. Furthermore, as 
d and A are assumed to be independent of the vertical coordinate, it is easy to check from the radial component of 
equation © that dBzfdz —>■ 0 (and hence —>■ 0). Therefore, the azimuthal and vertical components of equation 

©, at large R^, respectively lead to 
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dx 
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^ M-B,) = 0, (10) 

when the radial component of equation © turns out to be trivial. Subsequently the equation (|S]) reduces to 

i = 0- (11) 

Because of the choice of very large Rm (ideal MHD), there is a perfect flux freezing in the flow. Therefore, a steady 
advection of the vertical magnetic flux towards the center may lead to the decrease of (3, making it close to unity and 
further smaller, in a pure axisymmetric flow, even if the initial /3 was high. Hence, at some point, the back reaction 
of the field will inhibit accretion, depending on the geometry of the field lines. Although the physics of this process is 
not captured by the equations above and we also do not intend to discuss such physics, we will show below in §3 the 
effects of higher magnetic fields, in particular the inner edge (around the critical radius), in the entire flow structure. 
This is essentially to capture a situation after significant advection done with a certain field geometry. 

Henceforth, we will also neglect the second term in the parenthesis of equation ([5]). The set of equations o, @, 
©, 0, ©, ®, (nni) and m is essentially the modified version of the set of advective accretion disk equations in the 
presence of a large-scale magnetic effect, which are otherwise discussed in the literature in absence of it. Generally, in 
order to understand the hard state of accretion flows around black holes, the flow is assumed t o be purely hydrodynami c 
with the consideration of a turbulent viscosity arisen due to a weak magnetic field, i.e. MRI (iBalbus fc HawlewlllQQlll . 
Although, MRI is a largely accept ed idea, so far, in order to explain the origin of turbulen ce in accretion disks, there 
are some subtle issues with it (e.g. iMahaian fc Krishan II2008I : iMukhopadhvav et al. l[2Mlll including its applicability 
in colder disks. Therefore, transport of matter in disks is much more transparent through magnetic stresses, if the flows 
are embedded with a large scale field. Such magnetic stresses are considered here on the right hand side of the radial, 
azimuthal and vertical momentum balance equations. In a ddition, the magnetic heating due to abundant supply o f 
magnetic energy and the annihilation of the magnetic fields (|Bisnovatvi-Kogan fc Ruzmaikin 111974 IChoudhuri Ill998ll , 
which effects may be small however, is considered in the energy equation such that 


Q 


+ 

mag 


Q 


mag 


IQttx 


( 12 ) 


However, the other related terms in the momentum balance equations are neglected, again in comparison with the 
remaining terms in the respective equations, for the purpose of the present work. 

Therefore, even in the absence of turbulent viscosity (a = 0) and hence viscous stresses, magnetic stresses alone can 
help in transporting matter in the accretion flows. Such a consideration of large scale magnetic field and hence transport 
via m agnetic stress has, although been considered in circumstellar disks around young stars fsee lKdnigl fc Salmeroni 
I 2 OIII for a recent review), not yet been considered for advective accretion flows around black holes. 

Question may arise, if the magnetic field with plasma-/? > 1 is adequate enough to describe infall of matter in order 
to explain observation. We will show in the next section that the large-scale magnetic field, even with a significantly 
large plasma-/?, can describe advective accretion flows as efficiently as an a does. 


2.1. Solution procedure 

We have seven equations (excluding the vertical momentum balance equation, which assures the vertical magneto¬ 
static balance) and seven variables: d, A, P, p, B^, B^j,, B^, which we plan to solve along with the vertical magnetostatic 
balance condition. First, we plan to reduce d'd/dx in terms of other variables and the independent variable x alone 
(without any other derivatives), given by 
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As the advective accretion around black holes is necessarily transonic, the flow must pass through a critical radius 
where 
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when the variables with subscript ‘c’ indicate the respective values at that critical radius and the numerator of equation 
c© has to be zero for a continuous solution. At the critical radius, we also prescribe 


Bxc — Bye — Bze — \/ 47rp, 


—when Cs = sound speed ^ \ —, 

/aVS V P 


(15) 


so that the Alfven velocity is a fraction of sound speed therein. Although this is a simpler prescription, other choices 
do not change the picture, being addressed in this work, qualitatively. Note that a steady MHD flow would normally 
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have three critical points - fast m agnetosonic poin t, Alfven point and slow magnetosonic point - of which the Alfven 
point is not a true critical point (jGammie Ill999f l. The remaining two physically important distinct critical points, 
corresponding to fast and slow magnetosonic waves, collapse into a single point because of the assumptions made in 
equation m- We typically choose ^ 10 — 10^ in our various computations (see the figures). This is to capture a 
situation when magnetic pressure at the inner edge of the flow is not high enough to hinder radial infall of the matter, 
i.e. a situation with a weak back reaction of magnetic fields. Furthermore, from equation m , we can write 

XcBxc = constant = Cq = xB^ (16) 

which fixes the profile of B^ throughout the flow. 

These four conditions, along with the conditions that A = \k (when \k being the Keplerian angular momentum 
per unit mass) and d << 1 at the beginning of sub-Keplerian flow far away from the black hole, i.e. outer boundary, 
and ^ 1 at r ^ serve as important boundary conditions. Based on all the conditions, by solving the set of seven 
coupled differential equations, we can obtain the profiles for all variables including those of By and B^. Of course, 
then, one has to supply M, M, a, 7 and hence fms and fm for a flow. See lRaiesh fc Mukhonadhyav I (|2010aD for the 
solution procedure in details. 


3. SOLUTIONS 

Our main aim is to obtain the solutions of magnetized accretion flows. In other words, the aim is to understand 
how the large scale magnetic field (alone) can influence the mass transfer in accretion process, in particular in the the 
advective regime. This is important in the light of our ignorance of the origin of viscosity in accretion flows, which may 
be arisen from turbulence, as the molecular viscosity therein is well-known to be insignificant. Hence, our venture here 
is to investigate, if the large-scale magnetic field, however with large plasma-/?, can govern the same/similar transport 
of angular momentum as the well-known a-prescription does. Hence, we plan to understand the relative strengths 
between the magnetic stress tensors and the viscous stress tensors in order to control advective accretion flows. 

Our plan is to explore specifically two situations. (1) Flows with a relatively higher M and, hence, lower 7 , modelled 
around stellar mass black holes: such flows may or may not form Keplerian accretion disks. (2) Flows with a lower 
M and, hence, higher 7 , modelled around supermassive black holes: such flows are necessarily hot gas dominated 
advective (or advection dominated) accretion flows. 

3.1. Accretion around stellar mass black holes 

We choose M = 10Mq and M = 0.1 MEdd, when Mq and MEdd are the solar mass and the Eddington accretion 
rate respectively. However, this choice of M does not necessarily imply the flow to be purely Keplerian, rather ad- 
yective, which is iiideed, in general, under considerat ion. Such flows have temperature T > 10®K and p > 10“^gm/cc 
(|Sinha et al. 11200^ iRaiesh fc Mukhonadhyav Il2010a^ , which were extensively explore d in the context of the formation 
of shock in hot accretion disks and subsequent outflows and observed spectral states (IChakrabarti &: TitarchuiTIIlQQSl : 
IChakrabarti I[l99^ . A relation of 7 (and M) with the ratio of the gas and ra diation content of the flow and the cor¬ 
responding variations have been discussed bv iMukhonadhyav fc Dutta I (|2ni2l i. Therefore, following previous authors 
and for the convenience of comparison of the previous results without magnetic fields, we choose 7 = 1.335 along with 
the intermediate fms and fm- Note that while the results depend on the sign of S 2 , they do not depend on si and S 3 . 

Figure [T] compares accretion flows ( 1 ) in the presence of large scale magnetic field, but in the absence of a-viscosity: 
magnetic flow, and (2) in the presence of viscosity, but in the absence of large scale magnetic field: viscous flow. Here 
we consider to have the same direction as of A. It shows that the large scale magnetic field ^ 10"^ — lO^G, with its 
distribution in the inner edge of the flow defined by equation m, is adequately able to transport angular momentum, 
as viscous flows do with a = 0.017 and 0.012 respectively for nonrotating and rotating black holes. Figures [T^,b show 
that the disk sizes are the same for the respective viscous and magnetic sub-Keplerian flows with the above mentioned 
respective a-s. However, the transport of angular momentum takes place faster in magnetic flows, where the flows 
become quasi-spherical at larger radii than at the radii they become in the respective viscous flows. Away from the 
black hole, increases, which implies that the matter is prone to outflow from the outer region through the magnetic 
field lines extended outward direction. In a self-consistent model, including the flow variation in the vertical direction, 
the above features should have been appeared as the increasing magnetic field with the vertical coordinate. Such a 
model is planned to develop in future. In the present 1.5-dimensional magnetic flow model, when the magnetic stresses 
play the main role to remove angular momentum and hence to overcome the centrifugal barrier, as matter advances 
towards the black hole, the magnetic field decreases. Note that, as shown in Figs. [T}D,d, A and \B^\ both decrease 
towards the black hole: to overcome large A, the flow needs a large \B^\ and vice verse — they are the self-consistent 
solntions to each other. Note furthermore that the negative sign in in Fig. [TJd) is due to our choice of decreasing 
B^ with increasing z, i.e. negative S 2 , in the outer edge of the flow (which results in the same trend of the flow 
almost throughput, except very inner region). Figure[I](c) shows that our chosen regime of magnetic flows, allowing a 
steady infall of matter, corresponds to a relatively high plasma-/? (actually inverse of /? in shown). This furthermore 
renders lower magnetic pressures in respective flows compared to their maximum allowed value s based on the virial 
theorem/equipartition principle. As discussed in the previous works (e.g. IMukhonadhyav Il2003f) . a strong centrifugal 
effect is depicted in the Mach number profiles in Fig. [TJi (featured as slowing down the matter at around x = 20— 25) 
in the high angular momentum flows around nonrotating black holes, compared to the low angular momentum flows 
around rotating black holes. 
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Fig. 1. — (a) Mach number, (b) angular momentum per unit mass in GMjc, (c) inverse of plasma-/3, (d) azimuthal component of magnetic 
field in Gauss, when solid and dotted lines are for magnetic flows around Schwarzschild (a = 0, Ac = 3.2) and Kerr (a = 0.998, Ac = 1.8) 
black holes respectively, and dashed and long-dashed lines are for viscous flows around Schwarzschild (a = 0, ck = 0.017, Ac = 3.15) and 
Kerr {a = 0.998, o = 0.012, Ac = 1.8) black holes respectively. Other parameters are M = IOMq, M = O.lMEdd, 7 = 1-335,/^,is = fm = 
0.5, S2 = -0.5. 

Now we hypothesize that increases with the increasing z, i.e. positive S 2 , at the outer edge of the flow (which 
results in the same trend of the flow almost throughout, except very inner region). In this case, the right hand side of 
equation for a magnetic flow, completely flips sign compared to the negative S 2 case. Figure [2] shows that as the 
matter infalls, the toroidal component of magnetic field slows down the azimuthal motion of matter faster, making 
A = 0 and, subsequently, inverting the orientation of A. This is effectively due to the change in signs of magnetic 
stress components: and B^B^. Figure [2j) shows the variation of the magnitude of A, as the matter falls in. A 

is positive far away from the black hole, but it is negative close to the black hole. The location around A = 0 reveals 
a trough-like region in the flow. Hence, in the either sides of A = 0, there is a stronger centrifugal barrier which 
stores matter around A = 0 (due to the competition between radial and azimuthal flows). This region is prone to kick 
the matter out, producing outflows. Hence, if increases with z in the flow to start with, then as matter advances 
towards the black hole, a “potential well” forms to produce outflows. Note, however, that very close to a rotating black 
hole, matter will be under the influence of the black hole completely and hence A cannot have an opposite sign with 
respect to that of the black hole. Therefore, this solution is not valid very close to the rotating black hole. Indeed, the 
pseudo-Newtonian description is not applicable very close to the black hole. Nevertheless, the above solution implies 
a possibility of having such an origin of outflows in a magnetized accretion flow in the presence of a finite conductivity 
(when the field is not frozen with the matter, when the term associated with magnetic diffusivity in the induction 
equation is retained). Note that the plasma-/3 > 1 is maintained throughout the flows. 

This furthermore motivates us to check with such a possibility in viscous flows with viscosities a = 0.08 and 0.056 
respectively for nonrotating and rotating black holes. As shown in Figs. [lli,b, the angular momentum profiles in the 
respective magnetic and viscous flows appear to be similar, which furthermore makes the respective radial velocity 
profiles similar, unlike the previous cases, as shown in Fig. [I] with the positive A throughout. In the previous cases, 
the magnetic stresses are able to remove the angular momentum faster than the viscous stresses, in particular at a 
large distance from a Schwarzschild black hole, which is clearly understood from Figs. [TIi,b. Although the same is true 
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for a Kerr black hole, as the disk angular momentum itself is lower there, it does not effectively create any impact on 
the Mach number profiles. However, due to the choice of larger a, in the counterrotating cases, viscous stresses appear 
to be almost equivalent to the magnetic stresses and hence the radial velocity profiles in either of the respective flows 
appear similar. 

Let us now understand in more details, how the various components of magnetic stress are responsible for inflow 
and/or angular momentum transfer therein. Figures [3K,b show the variations of various components of the magnetic 
field as functions of the radial coordinate, around Schwarzschild and Kerr black holes, which are responsible for 
the various components of magnetic stress tensor. The profiles of field components and their magnitudes are partly 
dependent on their prescription given by equation (1151) . FigureH^ shows that the stress tensor component around 

a Schwarzschild black hole increases almost throughout as matter advances towards the black hole. This implies that 
the flow is prone to outflow through the field lines, which indirectly helps in removing the angular momentum, which 
furthermore renders its infall towards the black hole. However, very close to the black hole, B^B^ decreases, as indeed 
outflow is not possible in the near vicinity of the back hole, in particular, once the matter passes through the (inner) 
sonic point. By this radius, the flow angular momentum becomes very small which practically does not affect the 
infall. The magnitude of B^B^ decreases till the inner region of accretion flow, implying the matter to be spiralling 
out and hence removing the angular momentum. A larger \B^Bz\ at a larger radius implies a requirement of the 
removal of larger A therein. This automatically emerges from the self-consistent solutions of the set of equations. 
Nevertheless, close to the black hole, this effect reverses, rendering infall. Finally, the magnitude of B^B^ increases 
at a large and a small distances from the black hole (except around the transition radius), which helps infall, in 
the same way as the Shakura-Sunyaev viscous stress would do with the increase of matter pressure. However, at the 
intermediate zone, the transfer of angular momentum through B^B^ reverses and a part of the matter outflows. At the 
Keplerian to sub-Keplerian transition zone, flow/disk thickness increases, which effectively kicks the matter vertically, 
showing a decrease of B^B^. Most of these features remain unchanged for the flow around a rotating black hole, as 
shown in Fig. |3 }d. However, as a rotating black hole renders a stronger/efficient outflow/jet, here, except at the inner 
zone, BxB^ decreases throughout, which helps transferring the angular momentum inwards and kicking the matter 
outwards. Nevertheless, such a flow does not exhibit a high A either so that does not necessarily require an increasing 
BxB^ to remove A. Figures |4l:,d furthermore confirm that the above properties are invariant for the cases of B^ 
increased with increasing z. The only difference here is that B^B^j, and B^B^, have opposite signs with respect to the 
previous cases, when the components of magnetic field considered here are their respective averaged values. Note that 
all the components of magnetic stress tensor as functions of radius are determined by the associated components of 
magnetic field. The components of magnetic field are, however, determined by solving the underlying set of equations 
self-consistently with their prescription at the inner edge of the flow for the regime of interest. 

Figure El compares three cases of magnetic flows. (1) A counter rotating disk throughout, when B^ decreases with z 
almost throughout (solid line). (2) A disk having B^j, increasing with z almost throughout, which is corotating far away, 
but counterrotating close to the black hole (dotted line). (3) A disk having B^ increasing with z almost throughout, 
which is counterrotating far away, but corotating close to the black hole (dashed line). Note importantly in the latter 
two cases that the increasing B^ with disk height induces the change of handedness of the disk during the infall of 
matter. However, the profiles of Mach number and /3 practically appear similar in all the three cases. 

In Fig. 0 we compare the disk hydrodynamics between the magnetic flows with high and low magnetic fields. As 
expected, a flow with the higher magnetic field transports the angular momentum much faster, leading to a smaller 
sub-Keplerian flow. In other words, in the presence of a higher magnetic field, when the magnetic stresses are stronger, 
the Keplerian flow (when A = Xk) as well as the boundary between the Keplerian and sub-Keplerian flows are able 
to advance towards the black hole, shrinking the sub-Keplerian zone because of efficient angular momentum transfer. 
Figure[ 6 t shows that at a given radius the magnetic pressure, and hence the Alfven speed, is much larger in a flow with 
the higher magnetic field (but still /3 > 1 ). Naturally, such a high field magnetic flow would be equivalent to a viscous 
flow with much larger a, compared to the cases shown in Fig. [TJ An even higher magnetic field in the inner region 
would hinder any infall due to backreactions. Interestingly, in the radii close to the black hole, the sign of B^ becomes 
distinctly opposite than the outer region in the high magnetic field case. However, by this radius the required amount 
of angular momentum has already been transferred outwards in order to advance the matter close to the black hole 
and hence the change in sign of B^ does not create any physical impact onto the flow. If we vary the conditions chosen 
in equation (II5p . e.g., assume the components of magnetic field inequal, the qualitative picture remains unchanged 
— the magnetic stress in the presence of a large scale magnetic field could adequately transfer angular momentum. 
However, it is very important to note that if the strength of magnetic field and the corresponding value at the inner 
edge around the sonic radius would have been even higher, above a certain value, then the infall would no longer be 
possible. This is similar t o the situation when above a certa in value of A at the sonic radius in a given flow, the infall 
is no longer possible (see. lRaiesh fc Mukhonadhyav Il2010ah . 

If the flow is gas pressure dominated with larger 7 , then all the above basic features remain the same. Be it radiation 
or gas dominated, large scale magnetic stresses, yet ,5 > 1 , can transport angular momentum as efficiently as the a- 
prescription does. Nevertheless, below we discuss the effects of large scale magnetic field in the gas pressure dominated 
flows around a supermassive black hole. 

3.2. Gas dominated accretion around supermassive black holes 

The supermassive black hole at the centre of our galaxy, Sgr A*, presumably exhibits a gas pressure dominated, 
advection dominated, accretion flow. This motivates us to undertake this case, when we choose M = IO^Mq and 





log(x) 




log(x) log(x) 


Fig. 2.— Same as in Fig. ^ except S 2 = 0.5, when Ac = —3.2 and —1.65 for nonrotating and rotating magnetic flows respectively and 
a = 0.08, Ac = —3.08 and a = 0.056, Ac = —1.74 for nonrotating and rotating viscous flows respectively. 



Fig. 3.— Components of magnetic fleld: Bx (solid line), but normalized by 30 (dotted line), Bz (dashed line), for (a) Schwarzschild 
magnetic flow of Fig. ^ (b) Kerr magnetic flow of Fig. ^ 
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Fig. 4.— Components of magnetic stress: BxBfj, (solid line), BxBz (dotted line), Bfj^Bz (dashed line), for (a) Schwarzschild magnetic 
flow of Fig. ^ (b) Kerr magnetic flow of Fig. ^ (c) Schwarzschild magnetic flow of Fig. [2] (d) Kerr magnetic flow of Fig. [2] 

M = 10~‘^MEdd, and hence an appropriate 7 = 1.55. We furthermore choose fms = fm = 0.95, when strongly 
advective matter hardly has a chance to radiate photons. However, such a flow may also arise around a stellar mass 
black hole, e.g. microquasars. 

The basic features in Fig. [3 are similar as those for the stellar mass cases, as shown in Fig. [T] However, due to the 
gas dominance, and hence lower angular momenta, the Mach number profiles practically do not have any centrifugal 
barrier. Such flows are hotter, with T > 10^^K, and more quasi-spherical, compared to the radiation dominated flows, 
when a very small part of the dissipated heat can be radiated away. However, the most significant difference in these 
flows lies in their low magnetic fields, compared to those discussed in §3.1. This is due to the largeness of black hole 
masses in these flows, which leads to a much larger size of sub-Keplerian flows, when the dimensional flow size scales 
as M. As a result, due to the law of equipartition, the magnetic field decreases significantly compared to the cases of 
stellar mass black hole, as shown in Fig. [3i (as compared to those shown in Fig. [TJi). 

The magnetic flows around nonrotating and rotating, both the black holes, have their viscous counterparts with 
respective a = 0.011 and 0.0092. This furthermore confirms that even in gas dominated flows, the large scale magnetic 
field is able to transfer the angular momentum as efficiently as the a-prescription does with the most plausible values 
of a. 

Hypothesizing the increasing with increasing z at the outer edge of the sub-Keplerian flow, we obtain the same 
results as those for stellar mass black hole accretion flows described above, except at much lower magnetic fields. Like 
the stellar mass cases, here also a “potential well” forms which is featured in Fig. Hb, rendering the systems to have a 
zone for producing outflows. We do not repeat the detailed properties of it. Figures | 8 b,b also show the viscous flows 
resembling magnetic flows with a = 0.075 and 0.07 respectively for nonrotating and rotating black holes. 

3.3. Dependence on si, S 2 , S 3 

As defined in § 2 , si, S 2 , S 3 parametrize the scaling of the variations of i?^, B^ respectively in the vertical direction. 
This has be to considered because of our averaging the flows in the vertical direction, while the variation of magnetic 
field in the vertical direction has not been neglected, as has not been for P. Interestingly, the solutions practically do 
not depend on the choices of si and S3. However, with the increase of the magnitude of S2, which implies the increasing 
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Fig. 5.— Same as in Fig. [T] but all for magnetic flows around a black hole with a = —0.998, when the solid line is for S 2 = —0.5, Ac = 4.2 
and the dotted and dashed lines are for S 2 = 0.5 with Ac = 4.2 and —4.1 respectively. 

change of magnetic field with the vertical coordinate, the size of sub-Keplerian flow decreases. This is because, stronger 
the vertical variation of magnetic field, larger the change of the magnetic stresses in the vertical direction, on average 
faster the infall of matter is. This also argues for the faster rate of throwing the matter via outflows, when the 
outflows, in a more self-consistent 2.5-dimensional flow, are expected to plunge out via the magnetic field lines in 
the vertical direction. Hence, with the increase of magnetic field in the vertical direction, the system becomes more 
prone to outflow matter. Subsequently, a faster rate of outflow renders a faster removal of angular momentum and 
hence a faster rate of infall. As a result, the flow could remain Keplerian with the aid of adequate mechanisms of 
angular momentum transfer, till further inner region of the flow. Hence, the boundary between the Keplerian and 
sub-Keplerian flows advances towards the black hole. 

3.4. Interconnection between advection and magnetic field 

As the current J in the conducting fluid with conductivity a and electric field E is known to be J = cr (^E + v x , 

Faraday’s law of induction in the steady-state for axisymmetric accretion disks considered here, as given by equation 
©, can be recalled as 

V X ^-V X (ux h) =0. (17) 

For a flow with very large Rm (when Vm oc << 1), the z—component of equation (fTTll . averaged in z and integrated 
over (j), is given by 

J DBzxdcj) = constant = C, (18) 

when the constant can be identified as C = d/dt (J Bzdsx^) = d^/dt, where dsx(j> is the elementary surface area in 
the disk plane and $ is the magnetic flux. Therefore, equation (1181) fixes the relation between advection and B^, and 
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Fig. 6.— Same as in Fig. ^ but comparing between the flows with high (solid line) and low (dotted line) magnetic fields around a 
nonrotating black hole with S 2 = —0.5. 

hence the magnetic flux in the accretion flow. This also can be understood by recasting Faraday’s law of induction 
into 

( -* dA\ 

E + j = 0, when B = V x A, (19) 

which furthermore argues for 

f)A 

E = -VV - — -C, (20) 

when V is the Coulomb potential and (7 is a constant vector. Hence, for a steady axisymmetric accretion flow 

E^='^+xdB, = -C^. ( 21 ) 

a 

Therefore for Rm >> 1, is conserved. The constant or C can be fixed from a given boundary condition. In 

really, however, the flow is not expected to b e purely axisymmetri c and, hence, or C can also be mimicked as the 
contribution from non-axisymmetry. Earlier iLiibow et al. I (|1994h discussed a model of geometrically thin accretion 
flows in the presence of weak magnetic field, but assuming = 0 which is not true in general. 

The constraint on advection arisen in equation (|2T]) is clearly visible in the velocity profiles in Fig. [TJi with respect 
to the variations of the vertical component of magnetic field shown in Fig. |3l For the magnetic flow around a 
Schwarzschild black hole, d first increases steadily with the decrease of Bz at larger radii and subsequently matter 
tends to slow down due to centrifugal effect (with the relative decrease of infall rate dd/dx) with the increase of Bz 
until X = 10. Finally, matter plunges into the black hole steadily with a sharp decrease of Bz- For a rotating black 
hole, however, d steadily increases with the steady decrease of Bz almost throughout. Nevertheless, very close to the 
black hole, Bz slightly increases due to the spin effect of black hole, decreasing dd/dx slightly. This hints the power 
of black hole’s spin to plunge the matter out. 












12 



log(x) 





Fig. 7. — (a) Mach number, (b) angular momentum per unit mass in GMjc, (c) inverse of plasma-/!, (d) azimuthal component of magnetic 
field in Gauss, when solid and dotted lines are for magnetic flows around Schwarzschild (a = 0, Ac = 2.8) and Kerr (a = 0.998, Ac = 1.65) 
black holes respectively, and dashed and long-dashed lines are for viscous flows around Schwarzschild {a = 0, a = 0.011, Ac = 2.65) and 
Kerr (a = 0.998,0: = 0.0092, Ac = 1.52) black holes respectively. Other parameters are M = IO^Mq, M = 10~'^M^dd, 7 = 1-55, /t,is = 
fm = 0.95,82 = -0.5. 

4. DISCUSSION AND SUMMARY 

We have discussed the power of large scale magnetic field in advective accretion flows around black holes in oder to 
transport angular momentum, ena bling infall of matter. In a s impler Keplerian, self-similar model framework, such an 
investigation had been initiated bviBlandford &: Pavne I (|1982[) long ago, and in the cases of circumstellar disks around 
young stars (e.g. iKonigl fc Salmeron II2011I1 . such an approach has been explored. However, it remained unexplored, 
to the best of our knowledge, in the advective accretion disk around black holes, when it may exhibit a hard spectral 
state, until this work. Note that, often, only hard spectral states of disks are associated with the outflows/jets. 

We have found that the flows with plasma-/? > 1 exhibit adequate magnetic transport — as efficient as the a-viscosity 
with a = 0.08 would do. This is interesting as the origin of a (and the corresponding instability and turbulence) is 
itself not well understood. The maximum required large scale magnetic field is a few factor times lO^G in a disk 
around IQMq black holes and ^ lOG in a disk around W Mq supermassive black holes. The presence of such a field, 
in particular for a stellar mass black hole disk when the binary companion supplying mass is a Sun-like star with 
the magnetic field on average IG, may be understood, if the field is approximately frozen with the disk fluids (or the 
supplied fluids from the companion star remain approximately frozen with the magnetic field) or disk fluids exhibit 
large Reynolds number. Indeed, all the present computations are done at the limit of large Reynolds number, as really 
is the case in accretion flows, such that the term associated with the magnetic diffusivity in the induction equation can 
be neglected. The size of a disk around supermassive black holes is proportionately larger compared to that around 
a stellar mass black hole. Hence, from the equipartition theory, indeed the magnetic field is expected to be decreased 
here compared to that around stellar mass black holes. 

Is there any observational support for the existence of such a magnetic field, as required for the magnetic accretion 
flows discussed here? Interestingly, the polarizatio n measurements in th e hard state of Cyg X-1 imply that it should 
have at least lOmG field at the source of emi s sion dhaurent et al. Il201lll . In order to explain such high polarization, a 
jet model was suggested bv IZdziarski et alH (I20I4D . which requires a magnetic field ~ (5 — 10) x 10®G at the base of 
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Fig. 8.— Same as in Fig. [3 except S 2 = 0.5, when Ac = —2.8 and —1.62 for nonrotating and rotating magnetic flows respectively and 
a = 0.075, Ac = —2.65 and a = 0.07, Ac = —1.6 for nonrotating and rotating viscous flows respectively. 

jet and hence in the underlying accretion disk. Also, the magnetic field in the inner region of accretion disks around 
more than a dozen black holes has been found to be ver y high, based on a mod el relating the observed kinetic power 
of relativistic jet to the magnetic field of acretion disks ( Piotrovich et al. 1 IM1 . 

Different components of magnetic stress tensor have different roles: controls the infall in the disk plane, 

whereas B^B^ renders the flow to spiral outwards and, hence, outflow. Moreover, B^B^ helps to kick the matter out 
vertically. Larger the field strength, larger is the power of magnetic stresses. Interestingly, the magnitude of magnetic 
field decreases, as the steady-state matter advances towards the black hole. This is primarily because B,f,B^ (and also 
BxB^ for a rotating black hole) decreases inwards almost entirely in order to induce outflow. This furthermore reveals 
a decreasing as the output of self-consistent solutions of the coupled set of equations, which is reflected in the \B\ 
profile. 

In the present computations, we have assumed the flow to be vertically averaged without allowing any vertical 
component of the flow velocity (but keeping all the components of magnetic field). The most self-consistent approach, 
in order to understand vertical transport of matter through the magnetic effects which in turn leads to the radial 
infall of rest of the matter, is considering the flow to be moving in the ver tical direction from the disk plane as well. 
Such an attempt, in the absence of magnetic effects, was made earlier by iBhattacharva et al. I (j20intl in the model 
framework of coupled disk-outflow systems. In such a framework, the authors furthermore showed that the outflow 
power of the correlated disk-outflow systems increases with the increasing spin of black holes. Our future goal is now 
to combine that model with the model of present work, so that the coupled disk-outflow systems can be investigated 
more self-consistently and rigorously, when the magnetic field plays indispensable role in order to generate vertical 
flux in the three-dimensional flows. 
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